Lagrangian analysis of alignment dynamics for isentropic 
compressible magnetohydrodynamics 



J. D. Gibbon^ and D. D. Holmi'^ 
^ Department of Mathematics, Imperial College London SW7 2AZ, UK 

email: j.d.gibbon@ic.ac.uk and d.holm@ic.ac.uk 

^Computer and Computational Science, 

Los Alamos National Laboratory, 
MS D413 Los Alamos, NM 87545, USA 

email: dliolm@lanl.gov 

6 August 2006 



Abstract 

After a review of the isentropic compressible magnetohydrodynamics (ICMHD) equa- 
tions, a quaternionic framework for studying the alignment dynamics of a general fluid flow 
is explained and applied to the ICMHD equations. 
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1 Introduction to ICMHD 

Consistent with the topic of this special volume, this paper will deal with the equations of isen- 
tropic compressible magnetohydrodynamics (ICMHD) Not only are these sufficiently 
general to encompass most of the interest in this area but they also offer a significant challenge 
eyen without dissipation. The ICMHD equations determine the dynamics of a conducting fiuid 



1 



6 Aug 2006 



Lagrangian analysis of alignment dynamics for ICMHD 



2 



flow with magnetic field B satisfying divS = and current J = curlS. The equations for 
the fluid velocity field u, mass density p, specific entropy a, pressure p{p, a) and divergenceless 
magnetic field B are, respectively, the motion equation, Faraday's Law of frozen-in magnetic 
flux, isentropy along flow lines and the continuity equation : 

p^ = -Vp + JxB, ^ = curl (it X , (1.1) 

Do- Dp-^ . , , 

^ = 0, ^ = ,-div., (1.2) 

with the Lagrangian time derivative defined as 

D d , , 

The equation of state for specific internal energy e{p, a) and the Thermodynamic First Law 
yield 

so the specific internal energy changes only because of mechanical work. Combining Faraday's 
Law with the continuity equation yields 

DB 

-j^ - Bp-Vu = with Bp := p-^B . (1.5) 

This is the condition for the vector field Bp • V to be frozen into the fiow, i.e., 

D f „ d \ Dx 

-(b,.-)=0, along —■.= u. (1.6) 

We denote 

w:=p+\B'', (1.7) 

and set 

^:=B,.V, (1.8) 

with Bp defined in (|1.5() . In this notation, the ICMHD equations and H1.2|) transform using 
standard vector identities into 

Du dB 1 „ DB,, du 

^ = ^=p-Mivt.. (LIO) 

Dt Dt ^ ^ ' 

D/Dt and d/ds are defined in and (|1.8j) respectively. The first equation in expresses 

the motion equation in terms of the derivative along field lines, while the second expresses 

Faraday's Law H1.5() in terms of the derivatives D/Dt and d/ds. The invariance condition (|1.6|1 

for vector field d/ds = Bp ■ V along Lagrangian field lines for ICMHD implies equality of the 

following cross derivatives, 

D d d D 
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Of course, the invariant vector field d/ds = Bp ■ V may be applied to any fluid quantity. For 
example, applying d/ds to the equations in (|1.9|) and using the equahty of cross derivatives in 
d/ds and D/Dt yields the following exact nonlinear wave equations for ICMHD 

D\p~'B) d^B d( , ^ 

D'^u d'^u D / , \ , , 

-dW-1u^ = -mi" 

When linearized, these equations yield Alfven-sound waves. The present work emphasizes the 
alignment dynamics that is inherent in these equations, rather than their wave properties. 

1.1 Slow plus fast decomposition and nonlinear waves 

It seems natural to define slow and fast aspects of the ICMHD solutions in the sense of the 
Lagrangian time derivative. In particular, there are two relations for slow variables, 

^^0, and £1.0, (1.14) 

where /? is defined as 

P:=Bp-Va. (1.15) 

The latter follows from equality of cross derivatives, by computing 

D[j D ,^ ^ ^ D da d Da , ^ 

:=_(«,, V,) = __ = __=0. (1.16) 

Thus, the ICMHD equations preserve the projection of Bp on Vo" along Lagrangian flow lines, 
that is, along the path x{t) of a Lagrangian fluid particle determined from dx/dt = u{x, t). 
Equivalently, one finds the Eulerian conservation law, 

+ div(pPu) = , with pp = B-Va = dW(aB) , (1.17) 

ot 

upon combining (|1.16l) with the continuity equation for mass. 

Of course, the evolution of the component of Bp perpendicular to Vo" is not so simple! One may 
decompose Bp into components parallel and perpendicular to Vo" as 

iVcTpBp = /3VcT + 7 X Vcj = |Vo-p(SJI +S^) with 7 • Vcr = , (1.18) 

so that 

P = Bp-Va and 7 = Bp x Vcr . (1.19) 
One then computes the auxiliary equations. 



D _ _ . . Da 

and 



— -Vo- = - {Vuf ■ Va , that is, = - aju^^ (1.20) 



— \\/a\'- = -Va-S-Va (1.21) 
DV ^ ^ ' 

with the fluid strain-rate tensor defined as S = |(Vu + iS/u)^). Consequently, the evolution 

equation for the component of Bp perpendicular to Vu is determined from 

^ = -^(Bp X Va) = ^ X Va - Bp X [Vu^ ■ Va) (1.22) 
which is clearly neither fast nor slow. 
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1.2 Lagrangian dynamics of specific volume 

Two other likely candidates for fast variables are the specific volume and the velocity 
divergence divu, which satisfy 

-^=p-^diYU, (1.23) 

— divu = divF - iVuP . (1.24) 

Dt ' ' ^ ' 

The second of these equations is found using the identity 

D , /Du\ ,„ ,9 , , 

— divti = div - Vup 1.25 

Dt \Dt) ^ ^ ^ ' 

where |Vtip = u^^u'-. As a consequence of p.23() and (|1.24l) . the specific volume satisfies 
the "Lagrangian oscillation equation," 



P"' (^(diva)" - |VuI" + div(^— -/)-^Vti7jJ (1.26) 

upon substituting the definition of force F := ^ — p~^'Vw. Thus, the divergence of the force 
combines with compressibility of the flow to drive either Lagrangian oscillations or exponential 
variations of the specific volume, depending on the sign of the term on the right hand side of 
equation (|1.26j) . Of course, the sign of this term could be varying rapidly in some physical 
applications, and observations of such behavior would have corresponding implications for the 
dynamics of specific volume. In any case, the form of equation (|1.26|) is universal and one may 
conclude that the divergence of the total force drives the second Lagrangian time derivative of 
the specific volume. 

In what follows, we shall seek equations for alignment dynamics in MHD which have the same 
universal property as the Lagrangian oscillation equation (|1.26j) . In particular, we shall find 
that gradients of the total force (rather than its divergence) drive the Lagrangian dynamics of 
alignment in MHD of the frozen-in magnetic vector field Bp relative to its projection Bp ■ Vit 
onto the shear tensor Vit of the fluid velocity. 

1.3 Stretching and alignment 

The magnetic flux equation 

DBn du 

implies that the contravariant vector Bp undergoes stretching to the extent it aligns with the 
shear Vu. Likewise, this alignment evolves according to 

D'^Br, D du d Du d d / dB . \ , , 

/ = = = —F = —i p-^Vzu] , (1.28 

Dt^ Dt ds ds Dt ds ds\ ds ^ J ^ ' 

which recovers the nonlinear wave equation (|1.12|) above. 

We shall examine this equation from an alignment viewpoint, rather than as a wave propagation 
phenomenon. We begin by decomposing the vector Bp ■ Vu into its components parallel and 
perpendicular to Bp := Bp/\Bp\ 

Bp-Vu = aBp + xy< Bp, (1.29) 
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where a and x defined by 

a = Bp - (Bp-Vu), and x = Bp x {Bp ■ Vu) . (1.30) 

This decomposition is explained more fully in In fact, a is the Lagrangian amplification 

rate of the magnitude \Bp\ and x is the Lagrangian frequency of rotation of the unit vector Bp 
under the forcing by shear in the flux conservation equation ()1.27() . 

D\BJ , , , DBp 

-^ = a\Bp\ and —± = ^xBp. (L31) 

Note that no confusion should arise between this a and the one appearing in the a— dynamo 
equations! 

1.4 A mathematical framework for magnetic fluid alignment dynamics 

The quantity a defined in equation (|1.3()|) is the alignment of Bp with Vit, which determines 
the rates of change of magnitude \Bp\ whereas x is the misalignment, which determines the 
frequency of rotation of the direction Bp. The question asked in this situation is, "How long will 
Bp remain aligned with Vw, so it can continue to be stretched under the fluid shear?" Answering 
this question requires an analysis of Da/Dt and Dx/Dt using the alignment dynamics H1.28() . 

In principle, this analysis could be performed by direct computation and algebraic manipula- 
tion. However, the Lagrangian evolution equations for a and % happen to fit perfectly into a 
mathematical framework that was especially designed for analyzing orientation dynamics and for 
systematically interpreting the results. The key for recognizing this framework is to notice that 
the decomposition of a vector into its components parallel and perpendicular to another vector 
defines a type of product, or multiplication, that was first discovered by Hamilton [1]. This 
product reveals itself when we write the parallel-perpendicular vector decomposition equation 
H1.29() as though it were the pure vector components of a equation involving the four-component 
scalar-vector object (tetrad) [a, x] i^i the form, 

[0, Bp ■ Vu] = [0, a Bp + X X Bp] = [a, x] ® [0, Bp] with x-Bp = 0. (L32) 



This organization of the decomposition equation (|1.29j) summons the ® product defined by 

[pi, Qi] ® [p2, 92] = \p1p2 - qi ■ 92, piq2 + qiP2 + qi x 92] , (1.33) 

which is the multiplication rule that Hamilton invented for the field of quaternions A. The 
origin of this rule and its connection to the definition of a quaternion is given in N2.1I Of course, 
quaternions have everything to do with orientation [Sj IH) ■ The remainder of this paper sets up 
the quaternionic framework for alignment dynamics of a general fluid flow and applies it to the 
ICMHD equations. 

The plan of this paper is as follows: fJJ] summarizes the results of and shows how the evolution 
of the ortho- normal frame can be calculated in general terms. ^ then discusses the ICMHD 
equations (|1.9|) and in these terms, while fjl] examines alignment and growth properties 

to begin answering the question raised above for magnetic dynamos using these results. Finally, 
discusses the potential for other applications of this general method in compressible MHD 
turbulence. 
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2 Quaternions and Lagrangian alignment dynamics 
2.1 Background for quaternions 

The resurgence of practical interest in quaternions during the last two decades has been stim- 
ulated by progress in the computer animation and inertial navigation industries because of the 
ease with which quaternions handle moving objects undergoing three-axis rotations [SUnj. The 
tracking of aircraft and satellites and the animation of tumbling objects in computer graphics 
are obvious examples. Quaternion methods have also been recently been applied to the three- 
dimensional Euler equations for incompressible fluid motion j3 |Hj and to passive tracer particles 
transported by an underlying Lagrangian flow field; see [HI ITT] and references in The 
final result of these endeavors is that equations of motion can be derived for an ortho-normal 
frame. This 'quaternion frame' follows the evolution of particles in a Lagrangian flow whose 
evolution derives from the Eulerian equations of motion. 

Three-axis rotations lie at the heart of the definition of a quaternion ^3]. In terms of any 
scalar p and any 3-vector q, the quaternion c\ = [p, q] is defined (using Gothic fonts to denote 
quaternions) as 

3 

^ = \P, q] = Pl -"^Qif^i, (2-1) 

1=1 

where {ai, cr2, CT3} are the Pauli spin- matrices and / is the unit matrix. The relations between 
the Pauli matrices UiCTj = —SijI — eijkCk then give a non-commutative multiplication rule 



) q2 = \p1p2 - qi ■ 92, piq2 + qiP2 + qi x 92] • (2.2) 



It is easily demonstrated that quaternions are associative. In fact the individual elements of 
a unit quaternion provide the Cayley-Klein parameters of a rotation. This representation is 
a standard alternative to Euler angles in describing the orientation of rotating objects, as the 
books by Whittaker and Klein show jl4| II 5j . 

2.2 Quaternions and Lagrangian alignment dynamics in fluids 

A general quaternionic picture of the process of Lagrangian flow and acceleration in fluid dy- 
namics is explained in this section by considering the abstract Lagrangian flow equation 

Dw , , , , 

= a(x,t), (2.3) 

Dt ^ ^ V y 

whose Lagrangian acceleration equation is given in general by 

D'^w Da , , , , , 

b{x, t) . (2.4) 



Dt"^ Dt 

These are the rates of change of these vectors following the characteristics of the velocity gen- 
erating the path x{t) of a Lagrangian fluid particle determined from dx/dt = u{x, t). 
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{X2,t2) 



particle trajectory 




Figure 1: The dotted line represents the path of Lagrangian fluid particle (•) moving from {xi,ti) 
to {x2,t2)- The solid curves represent lines of constant w to which ttj is a unit tangent vector. The 
orientation of the quaternion-frame {w,x^, w x x^) is shown at the two space-time points; note that 
this is not the Frenet-frame corresponding to the particle path but to lines of constant w. 

Given the Lagrangian equation t\2.'A^ one defines the scalar Oa and the 3- vector Xa ^ 

aa = w'^{w-a), x„ = X a) , (2.5) 

in which w = ww with if = As observed in (|1.29() . the 3- vector a is decomposed into parts 
that are parallel and perpendicular to w as 



a = OaW X It) = [Oa, Xal ® [0> f^] 



(2.6) 



and thus the quaternionic product (|2.6|) is summoned in a natural manner. By definition, the 
growth rate Oa of the scalar magnitude w = \w\ obeys 



Dw 



a„w , 



while the unit tangent vector w = ww ^ satisfies 

Dw 



Dt 



X 



Now define two quaternions 



(\a = [aa, Xa\ and q;, = [ab, Xb\ 



(2.7) 



(2i 



(2S 



where ctb, Xb are defined as in (|2.5() for Ua, Xa with a replaced by b. Let ro = [0, w] be the pure 
quaternion satisfying the Lagrangian evolution equation 1)2. 3() with defined in (|2.9p . Then 
H2.3|) can automatically be re-written equivalently in the quaternion form 



'Dt 



[0, a] = [0, OaW X w] = qa ® tr . 



(2.10) 



Moreover, if a is Lagrangian-differentiable as in 1)2. 4|) then it is clear that a similar decomposition 
for b as that for a in H2.6|) gives 



D^tv 

Dt'^ 



[0, b] = [0, a^w X;, X w] = q;, ® tn . 



(2.11) 
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Using the associativity property, compatibility of (|2.11jl and ()2.1U|) implies that 

^ + qa®qa-qfej ®ra = 0, (2.12) 

which establishes a Riccati relation between the quaternions and q^ 

^ + qa®qa = qb. (2.13) 
From equation H2.13() there follows the main result of the paper: 

Theorem 1 The ortho-normal quaternion-frame F = {w, x^, wxx^) £ SO{3) has Lagrangian 
time derivatives expressed as 



Dw 
~Dt 

D{w X Xa) 
Dt 

DXa 

Dt 

where the Darboux vector T)a defined as 



'Da^w, (2.14) 

-Da^iw^Xa), (2.15) 
'DaXXa. (2.16) 



T>a = —W + Xa with Cb = Xb-{w X Xa) : (2.17) 

Xa 

is the angular frequency of rotation of the ortho-normal frame F. 
Remarks. 

• The frame orientation is controlled by the Darboux vector "Da = (cb/Xa, 0, Xa) which 
lies in the {w, x^) plane and is so named for its similarity to the Darboux vector in the 
Frenet-Serret equations for a space curve. Note that the vector X>a depends on Xb is 
independent of ab- 

• The frame dynamics equations 1)2. 14(1 - 1)2 .171) may also be re-written in matrix form by 
defining the 3x3 skew-symmetric matrix Ca G so(3) with entries [Ca]ij = — Cij^D^ as, 



(2.18) 



Proof : Finding an expression for the Lagrangian time derivatives of the components of the 
frame {w, w x x^, Xa)'^ requires the derivative of Xa- For this, one first recalls that the 3-vector 
b may be expressed in this ortho-normal frame as the linear combination 

w~'^b = abW + CbXa + dbi'^ >^ Xa) , (2.19) 
where Cb is defined in (|2.17j) and db = — {Xa ' Xb)- The 3-vector product Xb — w~^{w x b) yields 

Xb = Cb{w X Xa) - dbXa ■ (2.20) 






/ 





-Xa 





and Ca = 




Xa 





-Cb/x 




I 





Cb/Xa 
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To find the Lagrangian time derivative of Xa^ use tfie 3- vector part of tlie equation for the 
quaternion = [oa, Xa\ Theorem ^ 

^ = -2aaXa+Xb, ^ = -2aaXa-db, (2.21) 

where Xa — \Xa 

|. Using (IT^ and (IT^ there follows 

^ = CXaH^ X , R(^^!^_Xa)_ = ^^^-1^^ , (2.22) 

which gives equations 1)2 .141) - 1)2 .l?!) . ■ 

3 Application of the quaternionic alignment theorem to ICMHD 

The way is now clear to apply Theorem ^ of 321 to the ICMHD equations (|1.9jl and (|l.lUjl . First, 
we identify from (|r77|) and (fOH|) 

n, = B„ a = B,.Vu, ^ = ^^^^ = ' V (^) . (3.1) 

Although messy, this provides an explicit expression for b that depends upon the pressure p, 
Bp, J and their derivatives, as determined from the fluid evolution equations for u, p and a, 

h = BpV {-p'^Vp + JxBp) . (3.2) 

The alignment parameters {a, Xi (^bi Xb} identified as 

a = Bp-{Bp-Vu), X = BpX {Bp ■ Vu) , (3.3) 

ab = \Bp\-\Bp-b), Xb = \BprHBpXb). (3.4) 

These parameters appear in the vector decompositions, 

a = Bp-Vu = aBp + xx Bp, (3.5) 

and 

b = Bp-v{^^=a^Bp + Xb>^Bp. (3.6) 

The parameters Ub and Xb derive from Bp and from b in equation H3.2() at each time step. The 
vector b represents the coupling of the kinematic flow variables {Bp, Vu) to the gradients of the 
magnetic and thermodynamic forces. These identifications enter the two quaternions q = [a, x] 
and qb = [ab, Xb\ that satisfy the Riccati equation (|2.13)1 . The results of Theorem ^ for the 
Lagrangian time derivatives of the orthonormal frame F^ag = {Bp, Bp x x^x)'^ ^ 50(3) are 
then expressed as 

- VxBp, (3.7) 



Dt 

- = 'Dx{BpXx), (3.8) 



D{Bp X x) 

Wt - ^"^^^ 



Dx ^ . 
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where the Darboux angular velocity vector is defined as 

Cb 



X 



-Bp + x, with Cb = Xb- {Bp X x) , 



(3.10) 



which depends explicitly on the ICMHD force gradients through Xb is independent of ab- 

The frame dynamics equations ()3.7() - H3.10() for ICMHD may also be re-written in matrix form 
using the 3x3 skew-symmetric matrix C G so(3) with entries Cij = — eijkD^ as, 

/ 



DR 



mag 



Dt 



Remarks. 



CR 



mag 



where R, 



mag 




and C 





X 




-X 


Cb/X 




(3.11) 



DF, 



• The skew-symmetric matrix C = —^f^Rmag expresses the Lagrangian angular frequency 
of rotation of the magnetic orthonormal frame Rmag in terms of the magnitudes x ^^^^ ^b- 

• At a given moment in time, the flow lines of Bp may be constructed from its characteristic 
equations dx/ds = Bp{x). The Frenet-Serret (FS) equations describe how the orientation 
of an orthonormal frame changes along each flow line of Bp as a function of its shape pa- 
rameters, curvature and torsion. Applying equality of cross derivatives to equations ()3.11|) 
and (FS) yields a relation for the Lagrangian time derivatives of the shape parameters of 
a given flow line of Bp. 

Figure 1 now becomes 




{X2,t2) 



particle trajectory 




Figure 2: The solid curves represent magnetic field lines to which B is a unit tangent vector. The 
dotted line represents the path of a Lagrangian fluid particle (•) moving from {xi,ti) to (a;2,i2)- The 
orientation of the quaternion-frame (S, x, B x x) is shown at the two space-time points. 

4 Alignment and growth properties in ICMHD 

The growth rate a{x, t), defined in 1)3. 3|) . satisfies (see equations ()2.7|) and (12.81) ) 

^^^=a\Bp\. (4.1) 

Here a can take either sign and this is the key to how fast the magnitude \Bp\ increases (or 
decreases) at each point in the fiow. In contrast, the 3-vector x(a;, t) is the key to the alignment 
properties of the system, because it satisfies 

DBo - , , 

-^=XXB,^ (4.2) 
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As such, it can be interpreted as the swing rate of the unit vector Bp about a = Bp ■ Vu. 
Clearly, if Bp is aligned with a then X = and the quaternion q involves only the scalar a. 
Violent corkscrew-like motions of the magnetic field lines would be consistent with significant 
values of x s-nd such motions can therefore be regarded as a diagnostic for the misalignment of 
Bp with a. One of the messages of this paper is that putting a and x together as the quaternion 
q = [a, x\ is a natural way to approach this problem because the full quaternion q with X ^ 
is summoned whenever vortex or magnetic field lines bend or tangle. The Darboux vector is 
the angular frequency of rotation of the orthonormal frame F^ag = {Bp , Bp x x > x)"^) but is 
itself controlled by the gradients of forces in the expression ()3.2() for h. 

When written out in terms of a and x, the Riccati relation 1)2.13^ becomes 

Da 9 9 Dx , . . 

— = X^-a^ + ab, -j^ = -2ax + Xb ■ (4-3) 

Some years ago, these four equations were first expressed in this form for the incompressible 
Euler equations without any recourse to quaternions [T^. The quaternionic form first appeared 
in [S] (see [7|^] for a history). Equations 1)4. 3(1 appear to behave as Lagrangian ODEs driven 
by (]b = [o^bi Xb]- If the latter terms remain roughly constant equations (|4.3() can easily be shown 
to have two fixed points [12]; one has a negative value of a and the other has a positive value. 
The negative one is associated with an unstable spiral and the positive one with a stable spiral 
in the phase plane for this simplified ODE system. Equations (|4.3jl with constant q = [a^, Xb] 
have also recently been studied as a remarkably interesting kinematic model for the creation of 
non-Gaussian statistics in hydrodynamic turbulence .20^. This simple picture would seem not 
apply if the driving terms ab, Xb were to vary on the same time scales as a, X) or faster. This 
breakdown in applicability of equations ()4.3j) with constant q = [ab, Xb\ would be indicated if 
the force gradients in (|3.2() were observed to undergo rapid changes. 



5 Conclusion 

The quaternionic approach to Lagrangian frame dynamics developed here for ICMHD applies 
generally in fluid dynamics. Because the form of these equations is universal; that is, indepen- 
dent of the specific choice of forces, one may expect them to have many other applications. For 
example, in MHD one may use these equations to consider the effects on frame dynamics of 
rotation, or the effects of various kinds of subgrid-scale models, simply by identifying the corre- 
sponding expressions for the vectors a and b in equations (|2.3|) and (|2.4j) . Subgrid-scale models of 
MHD turbulence may be a particularly fruitful arena for these applications, especially for those 
models derived from Lagrangian averaging, because Faraday's Law for such models is preserved 
for the averaged field (21j. Depending as it does on Faraday's Law, the quaternionic method 
for MHD is fundamental, but it is also mainly kinematic. Thus its best role may be as a means 
of developing diagnostics for determining the effects of total force gradients. Therefore, addi- 
tional applications of the quaternionic approach may also be foreseen, as improved Lagrangian 
diagnostic methods are developed in the future. For example, future diagnostics may be able 
to distinguish between effects described by equations 1)4. 3(1 with fixed values of the quaternion 
c|b = [ctb, Xfe]) versus its self-consistent exact dynamics when the vector b is determined from the 
varying force gradients in equation (|3.2|) . 
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